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r^P ! Abstract 

This article introduces a metric for performance evaluation of medium access schemes in 
wireless ad hoc networks known as local capacity. Although deriving the end-to-end capacity 
of wireless ad hoc networks is a difficult problem, the local capacity framework allows us to 
quantify the average information rate received by a receiver node randomly located in the 
network. In this article, the basic network model and analytical tools are first discussed and 
applied to a simple network to derive the local capacity of various medium access schemes. Our 
goal is to identify the most optimal scheme and also to see how does it compare with more 
practical medium access schemes. We analyzed grid pattern schemes where simultaneous 
" transmitters are positioned in a regular grid pattern, ALOHA schemes where simultaneous 

transmitters are dispatched according to a uniform Poisson distribution and exclusion schemes 
where simultaneous transmitters are dispatched according to an exclusion rule such as node 
coloring and carrier sense schemes. Our analysis shows that local capacity is optimal when 
simultaneous transmitters are positioned in a grid pattern based on equilateral triangles and 
our results show that this optimal local capacity is at most double the local capacity of ALOHA 
based scheme. Our results also show that node coloring and carrier sense schemes approach 
the optimal local capacity by an almost negligible difference. At the end, we also discuss the 
^ . shortcomings in our model as well as future research directions. 
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1 Introduction 



This article presents the framework for deriving the local capacity in a wireless ad hoc network 
with various medium access schemes. The local capacity is defined as the average information rate 
received by a node randomly located in the network. The seminal work of Gupta & Kumar [T] 
and the later studies, e.g., [2H1] quantify the capacity in wireless networks in terms of scaling laws 
or bounds. These results are very important but they may not provide a deeper insight into the 
interrelationship between the performances of various medium access schemes in the network (e.g., 
the exact achievable capacity) and the possible network and protocol design issues (e.g., trade-offs 
involving protocol overhead versus performance of various medium access schemes) and network 
parameters. Apart from being of general interest, an additional advantage of local capacity, as 
compared to the results on scaling laws, is that it can be derived accurately. Therefore, the local 
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capacity framework can be used to get better insights into the designing of medium access schemes 
for wireless ad hoc networks. 

Medium access schemes in wireless ad hoc networks can be broadly classified into two main classes: 
continuous time access and slotted access. In this article, our main focus is on slotted medium 
access although many of our results can also be applied to continuous time medium access. Within 
slotted medium access category, we distinguish grid pattern, node coloring, carrier sense multiple 
access (CSMA) and slotted ALOHA schemes. In §2j we will briefly summarize our motivation and 
related works. In the following sections, we will concisely describe the analytical methods we have 
developed to derive the local capacity of the above mentioned medium access schemes in a wireless 
ad hoc network. Because these methods are developed independently, depending on the underlying 
design of each medium access scheme, ^3] will first give a brief overview of our network model and 
mathematical background of our parameters of interest. In sections HJ [5]and[6l we will discuss grid 
pattern, node coloring, CSMA and slotted ALOHA schemes respectively. First we will give a brief 
overview of the protocol overheads associated with each scheme and later, in each case, we will 
discuss the methods we used to evaluate local capacity of wireless ad hoc networks. In JTJ we will 
discuss the evaluation of these schemes and our results. The shortcomings and future extensions to 
our work and concluding remarks can be found in sections El and [9] respectively. 

2 Related Works 

In one of the first analyses on capacity of medium access schemes in wireless networks, [5] studied 
slotted ALOHA and despite using a very simple geometric propagation model, the result is similar 
to what can be obtained under a realistic SIR based interference model (non-fading, SIR threshold 
of 10.0 and attenuation coefficient of 4.0). Under a similar propagation model and assuming that 
all nodes are within range of each other, [6] evaluated CSMA scheme and compared it with slotted 
ALOHA in terms of throughput. [7J used simulations to analyze CSMA under a realistic SIR based 
interference model and compared it with ALOHA (slotted and un-slotted). For simulations, [7J 
assumed Poisson distributed transmitters with density A. Each transmitter sends packets to an 
assigned receiver located at a fixed distance of a\/A, for some a > 0. 

[8j|9] studied the transmission capacity, which is defined as the maximum number of successful 
transmissions per unit area at a specified outage probability, of ALOHA and code division multiple 
access (CDMA) medium access schemes. They assumed that simultaneous transmitters form a 
homogeneous Poisson point process (PPP) and used the same model for locations of receivers as 
in [7J. The fact that the receivers are not a part of the network (node distribution) model and are 
located at a fixed distance from the transmitters is a simplification. An accurate model of wireless 
networks should consider that the transmitters, transmit to receivers which are randomly located in 
their neighborhood or, in other words, their reception areas. Although, [5] analyzed the case where 
receivers are located at a random distance but it was assumed that a transmitter employs transmit 
power control such that signal power at its intended receiver is some fixed value. 

PPP only accurately models an ALOHA based scheme where transmitters are independently and 
uniformly distributed over the network area. However, under exclusion schemes, like node coloring 
or CSMA, modeling simultaneous transmitters as PPP leads to an inaccurate representation of 



2 



interference distribution. On the other hand, the correlation between the location of simultaneous 
transmitters makes it extremely difficult to develop a tractable analytical model and derive the 
closed-form expression for the probability distribution of interference. Some of the proposed ap- 
proaches are as follows. \±0\ (Chapter 18) used a Matern point process for CSMA based schemes 
whereas [TT] proposed to use Simple Sequential Inhibition (SSI) or an extension of SSI called SSIk 
point process. In [l2l[T3], interferers are modeled as PPP and outage probability is obtained by 
excluding or suppressing some of the interferers in the guard zone around a receiver. [H] analyzed 
transmission capacity in networks with general node distributions under a restrictive hypothesis 
that density of interferers is very low and, asymptotically, approaches zero. They derived bounds of 
high-SIR transmission capacity with ALOHA, using PPP, and CSMA using Matern point process. 

Other related works include the analysis of local (single-hop) throughput and capacity with slotted 
ALOHA in networks with random and deterministic node placement, and TDMA, in ID line- 
networks only [15] . [IE] determined the optimum transmission range under the assumption that 
interferers are distributed according to PPP whereas [17] gave a detailed analysis on the optimal 
probability of transmission for ALOHA which optimizes the product of simultaneously successful 
transmissions per unit of space by the average range of each transmission. [18] develops an ana- 
lytical model for analysis of saturation throughput of un-slotted CSMA with collision avoidance 
(CSMA/CA) in wireless networks. The numerical results of the analytical model are also verified 
using an event driven simulator. [19] uses queueing analysis to perform a comprehensive throughput 
and delay analysis of IEEE 802.11 MAC protocol. 

3 System Model 

We present our mathematical model in §3.11 and discuss the connection of local capacity with 
transport capacity in §3.21 

3.1 Mathematical Model and Assumptions 

We consider a wireless ad hoc network where an infinite number of nodes are uniformly distributed 
over an infinite 2D map. In slotted medium access, at any given slot, simultaneous transmitters in 
the network are distributed like a set of points 

S = {Zi, Z2, ■ ■ • , Z n , . . .} , 

where Zi is the location of transmitter i. The spatial distribution of simultaneous transmitters, i.e. 
the set S, depends on the medium access scheme employed by the nodes. Therefore, we do not 
adopt any universal model for the locations of simultaneous transmitters and only assume that, in 
all slots, the set S has a homogeneous density equal to A. 

We consider that each transmitter employs unit transmit power. The channel gain from node i to 
noise j is represented by 7^ such that the received power at node j is Pfiij, where is the transmit 
power of node %. We ignore multi-path fading or shadowing effects and assume that the channel 
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gain is solely determined by the distance and attenuation coefficient. Therefore, 



z j\ 



where |.| is the Euclidean norm of the vector and a > 2 is the attenuation coefficient. Therefore, 
the transmission from node i to node j is successful only if the following condition is satisfied 
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where N is the background noise (ambient/thermal) power and is the minimum signal to inter- 
ference ratio (SIR) threshold required for successfully receiving the packet. We assume that the 
background noise power, N , is negligible. Therefore, the SIR of transmitter i at any point z on 
the plane is given by 

\z — z-\~ a 
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We call the reception area of transmitter i, the area of the plane, A(z iy X, 0,a), where this trans- 
mitter is received with SIR at least equal to K. The area A(z iy X, 0, a) also contains the point z, t 
since here the SIR is infinite. The average size of A(z iy A, 0, a) is cr(A, 0, a), i.e., 

a{X,0,a) = E{\A{z h X,0,a)\) , 

where \A\ is the size of an area A. Note that cr(A, 0, a) does not depend on the location of transmitter 

Our principal performance metric is local capacity, hereafter referred to as capacity only, which 
is defined as the average information rate received by a receiver randomly located in the network. 
Consider a receiver at a random location z in the network and let N(z, 0, a) denote the number of 
reception areas it belongs to. Under general settings, following identity has been proved in [20] 



E(N(z,0,a)) = Xa(X,0,a) . (3) 

E(iV(z, 0, a)) represents the average number of transmitters from which a receiver, randomly located 
in the network, can receive with SIR at least equal to 0. Under the hypothesis that a node can 
only receive at most one packet at a time, e.g., when > 1, then N(z,0,a) < 1. The average 
information rate received by the receiver, c(z,0,a), is therefore equal to E(N(z, 0, a)) multiplied 
by the nominal capacity. Without loss of generality, we assume unit nominal capacity and we will 
derive 

c(z, 0, a) = E(N(z, 0, a)) = Xa(X, 0, a) . (4) 

We will derive the capacity in wireless ad hoc networks with grid pattern, node coloring, CSMA 
and ALOHA medium access schemes. We will also show that maximum capacity can be achieved 
with grid pattern schemes. Wireless networks of grid topologies are studied in, e.g., (2I1122] and 
compared to networks with randomly distributed nodes. In contrast to these works, we assume that 
only the simultaneous transmitters form a regular grid pattern. 
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3.2 Relationship of Local Capacity with Transport Capacity 



Gupta & Kumar pQ introduced the concept of transport capacity which measures the end-to-end 
sum throughput of the network multiplied by the end-to-end distance. A important aspect of their 
work and the following works on transport capacity is that it is not possible to compute the exact 
transport capacity in terms of network and system parameters and most of the results are in the 
form of bounds and scaling laws, i.e., the density of transport capacity scales as Cki\f\ bit-meters 
per second per unit area where C is the nominal capacity and k\ > depends on the medium access 
scheme and system parameters. We consider homogeneous traffic distribution in the network, i.e., 
all nodes have equal priority and generate traffic at the same rate. Therefore, if all nodes are capable 
of transmitting at C bits per second, the capacity of each node is Cki/y/X bit-meters per second. 

It is also shown in [1] that under general settings, the effective radius of transmission is fc 2 /vA for 
some k 2 > which also depends on the medium access scheme and system parameters. If each node 
transmits to a receiver which is randomly located within its effective radius of transmission or, in 
other words, its reception area, the information rate received by a receiver is constant and equal to 
Cki/k 2 bits per second. In this article, we evaluate the average of this information rate received by 
a receiver randomly located in the network, i.e., the local capacity. Note that, this capacity also 
incorporates the pre-constants associated with the scaling law, e.g. k\ and k 2 , and is independent of 
A as it is invariant for any homothetic transformation of the set of transmitters. Note that, we have 
abstracted out the multi-hop aspect of wireless ad hoc networks and this allows us to focus on the 
localized capacity. Our model also captures a realistic scenario of wireless ad hoc networks where 
each transmitter communicates with a receiver which is randomly located in its neighborhood or 
reception. 



4 Grid Pattern Based Schemes 

It can be argued that optimal capacity in wireless ad hoc networks can be achieved if simultaneous 
transmitters are positioned in a grid pattern. However, designing a medium access scheme, which 
ensures that simultaneous transmitters are positioned in a grid pattern, is very difficult because of 
the limitations introduced by wave propagation characteristics and actual node distribution. For 
this, location aware nodes may be useful but the specification of a distributed medium access scheme 
that would allow grid pattern transmissions is beyond the scope of this article. 

In this section, first we will present the analysis we used to investigate the optimality of grid 
patterns based medium access schemes. Later, we will also discuss the analytical method we used 
to to analyze the capacity of such schemes. Grid pattern based medium access schemes may have 
no practical implementation but their evaluation is interesting in order to establish an upper bound 
on the optimal capacity in wireless ad hoc networks. In the later sections, we will also discuss more 
practical medium access schemes. 
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4.1 Optimality of Grid Pattern Based Schemes 



In this section also, we consider that an infinite number of transmitters are uniformly distributed 
like a set of points 

S = {Zi, Z2, ■ ■ ■ , Z n , . . .} , 

on an infinite 2D plane. The location of transmitter i is denoted by Zi and the center of the plane 
is at (0,0). 

In order to simplify our analysis, we define a function as 



Si{z) = 
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where a > 2. The function Si(z) is similar to the SIR function Si(z), in (jSJ), except that the 
summation in the denominator factor also includes the numerator factor. In order to simplify the 
notations, we will remove the reference to z when no ambiguity is possible. 

We also define a function f(si) which can be continuous or integrable. For instance, here we will 
use 



for some given Note that, in this case, the function is not continuous but we will not bother 
with this. In the following discussion, we can consider without loss of generality that the value of 
0' is related to f3 by 




Therefore, if transmitter % is received successfully at location z (or in other words, with SIR at least 
equal to 0, i.e., Si(z) > (3), then Si(z) > f3' and /(sj) is equal to 1. 

We also assume a virtual disk on the plane centered at (0, 0) and of radius R. This allows us to 
express the density of set S, v(S), in terms of the number of transmitters covered by the disk of 
radius R or area nR 2 , when R approaches infinity, and it is given by a limit as 



We denote h(z) = J2i f( s i)- Note that, h(z) is equal to the number of transmitters which can be 
successfully received at z and its maximum value shall be 1 if /3 > 1. 

We define E(/i(z)) by the limit 



Note that, the integration is over an infinite plane or, in other words, over the disk of radius R 
where R approaches infinity. Also note that the notations are simplified by taking dxdy equal to 
dz 2 . We denote the reception area of an arbitrary transmitter i as 



f(Si) = l Si {z)>P' , 
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and we have 

E(h(z)) = l™ iE 1 !^^ = K$)E(*i) 



R->oo irR 2 

with 



E(o-j) = lim - V" o-j . 



n— >oo -n, 

i<n 

As i? approaches infinity, n, i.e., the number of transmitters in the set S, covered by the disk of 
radius R, also approaches infinity. 

Our objective is to optimize ~E(h(z)) whose definition is equivalent to the definition of E(N(z, K, a)) 
and therefore capacity, c(z,K,a), as well in expressions (j3j) and (Tj0) respectively 

4.1.1 First Order Differentiation 

We denote the operator of differentiation w.r.t. z% by V«. For % ^ j, we have 

Z - Zi 



and 



Therefore, 



\z — zA 2 



T7 I 2 ^ Z Z, 

V,s,- = a(s 4 — s 



l z - & H 
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Although, we know that 

J h(z)dz 2 = oo , 

we nevertheless have a finite V$ / h(z)dz 2 . In other words, the sum ^ • VjCij converges for all i. 
Lemma 1. For all j in S, 

i 

Indeed this would be the differentiation of Oj "when all points in S are translated by the same vector. 
Similarly, 

J2 Vi J h (z)dz 2 =0 . 
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Theorem 1. If the points in the set S are arranged in a grid pattern then 



Vj / h(z)dz 2 = ^ViO-j = , 



and grids patterns are locally optimal. 



Proof. If S is a set of points arranged in a grid pattern, then 



V, / h(z)dz 2 = J^ViO-j , 



would be identical for all % and, therefore, would be null since 



Vj / h(z)dz 2 = . 



We could erroneously conclude that, 

- all grid sets are optimal and 

- all grid sets give the same E(h(z)). 

In fact this is wrong. We could also conclude that E(<7j) does not vary but this will contradict that 
u(S) must vary. The reason of this error is that a grid set cannot be modified into another grid set 
with a uniformly bounded transformation, unless the two grid sets are just simply translated by a 
simple vector. □ 

However, we have proved that the grid sets are locally optimal within sets that can be uniformly 
transformed between each other. In order to cope with uniform transformation and to be able to 
transform a grid set to another grid set, we will introduce the linear group transformation. 

4.1.2 Linear Group Transformation 

Here, we assume that the points in the plane are modified according to a continuous linear transform 
M(t) where M(t) is a matrix with M(0) = I, e.g., M(t) = I + tA where A is a matrix. 

Without loss of generality, we only consider cq, i.e., the reception area of the transmitter at z 
which can be located anywhere on the plane. Under these assumptions, we have 





= B, 
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such that 



D 



DyX Dyy 



Therefore, we can write the following identity 

d 



tr ( A M a »)4 ffo(t ' A) 

where a (t, A) is the transformation of a under M(t), i.e., 

a (t, A) = det(I + tA)o- . 
We assume that M(t) = (1 + t)I with A = I, i.e., the linear transform is homothetic. 
Theorem 2. D is symmetric and tr(D) = 2a . 

Proof. Under the given transform, 

<7 (t,A) =<7 (M) = (l + t) 2 a . 

As a first property, we have 

tr(D) = 2a , 

since the derivative of <To w.r.t. identity matrix I is exactly 2o"o (i.e., tr(A T D) = tr(D) = <Tq(0,I) = 2<To). 
The second property that D is a symmetric matrix is not obvious. The easiest proof of this property 
is to consider the derivative of o" w.r.t. the rotation matrix J given by 



t=o 



J = 



-1 

1 



which is zero since J is the initial derivative for a rotation and reception area is invariant by rotation. 
Therefore, 

tr(J T D) = D yx — D xy — , 
which implies that D is symmetric. □ 

Note that D can also be written in the following form 

D = ^2 z i ® Vitro = J dz 2 ^Zi® Vif{s Q ) . 



Let T be defined as 



such that 



T = J dz 2 ^y z - ® Vi/(s ) , 



D= f ^^®V i /(s )^ 2 -T . 
J i 

The purpose of these definitions will become evident from theorems 3 and 4. 
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Theorem 3. We will show that JJ2 { z <8> Vif(s )dz 2 is equal to a I and, therefore, 

D = a I- T . 

We will also prove that T is symmetric. 

Proof. From the definition of T, we can see that the sum Y1A Z ~ z i) ® ^if( s o) leads to a symmetric 
matrix since 

T = a f /'( So ) (A—^(z ~ z ) ®(z- zo) + 
J K \z — ZoY 

S ° Sl , 2 {z - Zi) <S>(z- ZiYjdz 2 , 



E 



Z - Zi 


















and the left hand side is made of (z — Zi) <g> (z — Zi) which are symmetric matrices. This implies 
that T is also symmetric. 

We can see that 

E V */( s o) = -W(so) , 

i 

and using integration by parts we have 

/ J2 Z ® Vif(s )dz 2 = -lx 
J i 

J x£f(s )dxdy J x-§-f(s )dxdy 
Jy-^f(s )dxdy J y^f(s )dxdy 

which is symmetric and equal to <ToI. The sum/ difference of symmetric matrices is also a symmetric 
matrix and, therefore, D is a symmetric matrix and D = ool — T. □ 

Now, we will only consider grid patterns and, by virtue of a grid pattern, we can have 

Efo) = o- = J f(s )dz 2 , 

and E(h(z)) = v(S)a . Under homothetic transformation, v(S) and a are transformed but z/(5)cr 
remains invariant. 

Theorem 4. If the pattern of the points in set S is optimal w.r.t. linear transformation of the set, 
D = a I and T = 0. 

Proof. The derivative of a w.r.t. matrix A is exactly equal to D. Similarly, under the same 
transformation 

Fi U(S) = det(I + At) u(S) ' 

and for A = I, it can be written as 



u'(S)(t,l) = u(S)/(l + t) 2 . 
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Figure 1: Square, Hexagonal and Triangular grid patterns. The arrows (blue and red) represent 
the invariance of Eigen values w.r.t. isometric symmetries of the grids. 
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Figure 2: Computation of the reception area of transmitter i. 



In any case, the derivative of v(S) w.r.t. matrix A is exactly equal to —lu(S). We also know that 
if the pattern is optimal w.r.t. linear transformation, the derivative of u(S)<yo w.r.t. to matrix A 
shall be null. This implies that 

i/(S)D - lu(S)a = , 



which leads to D = ctqI and T = 0. 



□ 



We know that T is symmetric and T = 0. Thus, tr(T) = 0, i.e., Eigen values are invariant by 
rotation. When a grid is optimal, we must have T = 0. In any case, the matrix T must be invariant 
w.r.t. isometric symmetries of the grid. On 2D plane, the grid patterns which satisfy this condition 
are square, hexagonal and triangular grids. The square grid is symmetric w.r.t. any horizontal or 
vertical axes of the grid and, in particular, with rotation of 7r/2 represented by J. Therefore, the 
Eigen system must be invariant by rotation of tt/2. This implies that the Eigen values are the same 
and therefore null since tr(T) = 0. Same argument also applies for the hexagonal grid with the 
invariance for 7r/3 rotation and for the triangular pattern with invariance for 27r/3 rotation. 



4.2 Reception Areas 

The simultaneous transmitters, i.e., the set S is a set of points arranged in a grid pattern. We 
consider that, for every slot, the grid pattern is the same modulo a translation. We have covered 
grid layouts of square, hexagonal and triangle as shown in Fig. [TJ Grids are constructed from d 
which defines the minimum distance between neighboring transmitters and can be derived from the 
hop-distance parameter of a typical TDMA-based scheme. The density of grid points, A, depends 
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on d. However, the capacity, c(z, (3, a), is independent of the value of d or, for that matter, A as it 
is invariant for any homothetic transformation of the set of transmitters. 

Our aim is to compute the size of the reception area, A(Zi, A, (3, a), around each transmitter i. By 
consequence of the regular grid pattern, all reception areas are the same modulo a translation (and 
a rotation for the hexagonal pattern), and their surface area size, <r(A, /3, a), is the same. 

If C(zi,/3,a) is the closed curve that forms the boundary of A(zi, X, (3,a) and z is a point on 

C(zi, (3, a), we have 

a(X,/3,a) = ^ J det(z - z h dz) , (5) 

where det(a, b) is the determinant of vectors a and b and dz is the vector tangent to C(zi,f3,a) at 
point z. det(z — Zi,dz) is the cross product of vectors (z — Zj) and dz and gives the area of the 
parallelogram formed by these two vectors. Note that, (jSJ) remains true if Zi is replaced by any 
interior point of A(zi, X, (3,a). 

The SIR Si(z) of transmitter i at point z is given by (TSJ). We assume that at point z, Si(z) = (3. 
On point z we can also define the gradient of Si(z), VSi(z), 



SAz) 
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Note that VSi(z) is inward normal to the curve C(zi, (3, a) and points towards z%. The vector dz is 
co-linear with J .^^j, where J is the ant i- clockwise rotation of 3n/2 (or clockwise rotation of 7r/2) 
given by 

1 
-1 



J 



Therefore, we can fix dz = J m^f^ fit anc i i n © 



det(z — Zi, dz) 



. VSjjz) xi 
z - Zi) x J , :, 5t 



-(Z - Zi 



\VSi{z 
VS t (z) 
\VSi(z) 
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where St is assumed to be infinitesimally small. The sequence of points z(k) computed as 

z(0) = z 

/, -\ \ /, x t ySAzik)) . 
z fc + 1 ) = z (k + J ' ft, 
|V5i(2;(fc))| 

gives a discretized and numerically convergent parametric representation of C(zi, (3, a) by finite 
elements. 



Therefore, (El) reduces to 



a(X, (3, a) 



- Zi 



VSj(z(k)) 
■\VSi(z(k))\ 



St 



(6) 
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assuming that we stop the sequence z(k) when it loops back on or close to the point z. 

The point, z(0) = z, can be found using Newton's method. First approximate value of z, required 
by Newton's method, can be computed assuming only one interferer nearest to the transmitter % as 
discussed in the Appendix. The negative sign in (Q is automatically negated by the dot product 
of vectors (z(k) — z^) and VSi(z(k)). 

4.3 Capacity 

The expression to derive the capacity is 

c(z,P,ct) = E(N(z,P,a)) = N(z,P,a) = Xa(X,P,a) , 
where a(X, (3, a) is computed using the above described method. 

5 Medium Access Schemes Based on Exclusion Rules 

As we mentioned earlier, an accurate model of interference distribution, in case of exclusion schemes, 
is very hard to derive because of the correlation in the locations of simultaneous transmitters, 
e.g., simultaneous transmitters separated by a certain minimum distance. We also mentioned a 
few approaches to model distribution of transmitters in exclusion schemes including Matern point 
process PUCE] or, more recently, SSI point process [TT] . 

Classical Matern point process is a location dependent thinning of PPP such that the remaining 
points are at least at a certain distance b > from each other. Instead of distance, [TUHH] used 
received power level, a function of distance, in order to model CSMA based schemes. However, [TT] 
showed that Matern point process may lead to an underestimation of the density of simultaneous 
transmitters and proposed an SSI point process in which a node can transmit or, in other words, 
can be added to the set of simultaneous transmitters S if it is at least at a certain distance b > 
from all active transmitters (i.e., the transmitters already in the set S). SSI point process can 
be used to study node coloring schemes but may result in a flawed representation of CSMA based 
schemes. In case of CSMA based schemes, a node senses the medium to detect if the signal level 
of interference is below a certain threshold and only transmits if it remains below that threshold 
during the randomly selected back-off period. Therefore, the decision of transmission depends on all 
nodes which are already active and transmitting on the medium. In order to address this inaccuracy 
in SSI point process model, [TT] proposed an SSIk point process which ensures that a node, before 
transmitting on the medium, takes into account the interference from k nearest active transmitters. 
However, very few analytical results are available on SSI and SSI^ point processes and most of the 
results are obtained via simulations. Therefore, we will use Monte Carlo simulation along with the 
analytical method proposed in §4.21 to compute the capacity of node coloring and CSMA based 
schemes. 

In this section, we will discuss the models of node coloring and CSMA based schemes which we will 
employ in our Monte Carlo simulation later. In the following discussion, the set of all nodes in the 
network is N. In practical implementation, this set is finite but in theory, it can be infinite but with 
a uniform density. 
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5.1 Node Coloring Based Schemes 

Node coloring schemes use a managed transmission scheme based on time division multiple access 
(TDMA) approach. The aim is to minimize the interference between transmissions that cause 
packet loss. These schemes assign colors to nodes that correspond to periodic slots, i.e., nodes that 
satisfy a spatial condition, either based on physical distance or distance in terms of number of hops, 
will be assigned different colors. This condition is usually derived from the interference models of 
wireless networks such as unit disk graph (UDG) models. For example, in order to avoid collisions 
at receivers, all nodes within k hops are assigned unique colors. Typical value of A; is 2. A few 
practical implementations of node coloring schemes are as follows. In [23J, authors proposed coloring 
based on RAND, MNF and PMNF algorithms. In RAND, nodes are colored in a random order 
whereas MNF and PMNF prioritize nodes according to the number of their neighbors. NAMA [21] 
colors the nodes, in 2-hop neighborhood, randomly using a hash number. In SEED EX (25]) nodes 
use random seed number in 2-hop neighborhood to randomly elect the transmitter. FPRP [26] is 
a five-phase protocol where nodes contend to allocate slots in 2-hop neighborhood. DRAND is the 
distributed version of RAND and article [27] shows its better performance as compared to SEEDEX 
and FPRP. [28J proposes a joint TDMA scheduling/load balancing algorithm for wireless multi-hop 
networks and shows that it can improve the performance in terms of multi-hop throughput and 
fairness. Most of these schemes use unit disk graph based interference model. However, success of 
a transmission depends on whether the SIR at the receiver is above a certain threshold. [29] is an 
example of a node coloring scheme which uses SIR based interference model. Note that, extremely 
managed transmission scheduling in node coloring schemes has significant overhead, e.g., because 
of the control traffic or message passing required to achieve the distributed algorithms that resolve 
color assignment conflicts. 

In this article, instead of considering any particular scheme, we will present a model which ensures 
that transmitters use an exclusion distance in order to avoid the use of the same slot within a 
certain distance. This exclusion distance is defined in terms of euclidean distance d which may be 
derived from the distance parameter of a typical TDMA-based scheme. Therefore, a slot cannot 
be shared within a distance of d or, in other words, nodes transmitting in the same slot shall be 
located at a distance greater or equal to d from each other. 

Following is a model of node coloring schemes which constructs the set of simultaneous transmitters, 
S, in each slot (this is supposed to be done off-line so that transmission patterns periodically recur 
in each slot). 

1. Initialize M = M and 5 = 0. 

2. Randomly select a node Sj from M. and add it to the set S, i.e., S = S U {s^}. Remove Sj 
from the set A4. 

3. Remove all nodes from the set M. which are at distance less than d from Sj. 

4. If set Ai is non-empty, repeat from step 2. 

Above described steps model a centralized or distributed node coloring scheme which randomly 
selects the nodes for coloring while satisfying the constraints of euclidean distance. Under given 
constraints, this model also tries to maximize the number of simultaneous transmitters in each slot 
and should give the maximum capacity achievable with any node coloring scheme which may not 
prioritize the nodes for coloring, e.g., [27J. 
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5.2 CSMA Based Schemes 



As compared to managed transmission schemes like node coloring, CSMA based schemes are simpler 
but are more demanding on the physical layer. Before transmitting on the channel, a node verifies 
if the medium is idle by sensing the signal level. If the detected signal level is below a certain 
threshold, medium is assumed idle and the node transmits its packet. Otherwise, it may invoke a 
random back-off mechanism and wait before attempting a retransmission. CSMA/CD (CSMA with 
collision detection) and CSMA/CA (CSMA with collision avoidance), which is also used in IEEE 
802.11, are the modifications of CSMA for performance improvement. 

We will adopt a model of CSMA based scheme where nodes contend to access medium at the 
beginning of each slot. In other words, nodes transmit only after detecting that medium is idle. We 
assume that nodes defer their transmission by a tiny back-off time, from the beginning of a slot, 
and abort their transmission if they detect that medium is not idle. We also suppose that detection 
time and receive to transmit transition times are negligible and, in order to avoid collisions, nodes 
use randomly selected (but different) back-off times. Therefore, the main effect of back-off times is 
in the production of a random order of the nodes in competition. 

For the evaluation of the performance of CSMA based scheme, we will use the following simplified 
construction of the set of simultaneous transmitters S. 

1. Initialize M = Af and <S = 0. 

2. Randomly select a node s« from M. and add it to the set S, i.e., S = S U {sj}. Remove Sj 
from the set A4. 

3. Remove all nodes from the set A4 which can detect a combined interference signal of power 
higher than 9 (carrier sense threshold), from all transmitters in the set S, i.e., if 

J2\zi-z 3 \- a >e , 

remove s,- from A4. Here, Z; t is the position of Sj and | Z^ Z j I IS the euclidean distance between 
Si and Sj. 

4. If set A4 is non-empty, repeat from step 2. 

These steps model a CSMA based scheme which requires that transmitters do not detect an inter- 
ference of signal level equal to or higher than 9, during their back-off periods, before transmitting on 
the medium. At the end of the construction of set S, some transmitters may experience interference 
of signal level higher than 9. However, this behavior is in compliance with a realistic CSMA based 
scheme where nodes, which started their transmissions, or, in other words, are already added to the 
set S do not consider the increase in signal level of interference resulting from later transmitters. 

5.3 Reception Areas 

We are not aware of an analytical closed-form expression for the probability distribution of signal 
level with node coloring or CSMA based schemes. Consequently, we do not have a closed form 
expression for the average size of the reception area of an arbitrary transmitter with these schemes. 
Therefore, it is evaluated via Monte Carlo simulation using the analytical method of §4.21 The 
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value of d, in case of node coloring scheme, or 9, in case of CSMA based scheme, can be tuned to 
obtain an average transmitter density of A. 

5.4 Capacity 

Similarly 

c(z, (3, a) = E(N(z, (3, a)) = Acr(A, 0, a) , 

is also computed via Monte Carlo simulation. The capacity, c(z,/3,a), is invariant for any homoth- 
etic transformation of A and, therefore, it is also independent of the values of medium access scheme 
parameters 9 or d. 

6 Slotted ALOHA Scheme 

In slotted ALOHA scheme, nodes do not use any complicated managed transmission scheduling 
and transmit their packets independently (with a certain medium access probability), i.e., in each 
slot, each node decides independently whether to transmit or otherwise remain silent. Therefore, 
the set of simultaneous transmitters, in each slot, can be given by a uniform Poisson distribution 
of mean equal to A transmitters per unit square area [911171120] . 

We can write ([T]) as 

\z-Zi\~ a > $Yl\ z ~ z i\~ a 
or W{z, {zi}) > pW{z,S- {Zi}), where W{z,S) = £ |* - Zj\~ a . 

6.1 Distribution of Signal Levels 

Figure [3] shows the function W(z,S) for z varying in the plane with S an arbitrary set of Poisson 
distributed transmitters. Figure [3] uses a = 2.5. It is clear that closer the receiver is to the trans- 
mitter, larger is the SIR. For each value of (3 we can draw an area, around each transmitter, where 
its signal can be received with SIR greater or equal to (3. Figure H] shows reception areas for the 
same set S, as in Fig. EJ for various values of (3. As can be seen, the reception areas do not overlap 
for (3 > 1 since there is only one dominant signal. For each value of [3 we can draw, around each 
transmitter, the area where its signal is received with SIR greater or equal to (3. Our aim is to find 
the average size of this area and how it is a function of A, f3 and a. 

W(z,S) depends on S and hence is also a random variable. The random variable W(z, S) has a 
distribution which is invariant by translation and therefore does not depend on z. Therefore, we 
denote W(A) = W(z, A). Let w(S) be its density function. The set S is given by a 2D Poisson 
process with intensity A transmitters per slot per unit square area and Laplace transform of w(S), 
w(9, A), can be computed exactly. The Laplace transform, w(9, A) = exp( J (e~ 9r " — l)rdr), satisfies 
the identity 

w(9, A) = exp ( -7rAr(l - 




16 



10"" o 



Figure 3: Signal Levels (in dB's) for a random network with attenuation coefficient a = 2.5. 




Figure 4: Distribution of reception areas for various value of SIR. (3 = 1,4, 10 for situation of Fig. 
® 



where T(.) is the Gamma function. 

Note that, in all cases, w(9, A) is of the form exp(— AC# 7 ) where 7 = -, and the expression of C in 
case of 2D map is C = 7iT(l — 7). 

From the above formula and by applying the inverse Laplace transformation, we get 



P(W(A) < x) 



+ioo 

1 fw(9,X) 



e dx d6 



2m J 9 

-ioo 
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Expanding w(6, A) = £ ^ ^^-^ r^7 ) we get 

n>0 



+ioo 

P(W(A) <X) = ±_yt£^l I Qni-l e 9x d9 

v v ' ' 2nx^ n\ J 



By bending the integration path towards negative axis 



CO 



— / e^e^de = / e n ~<- l e- ex d6 

2m / 2in 



o 



7T 

we get, 



T(n7)x 



PWI < l) = ^^ir Wl -.. (7) 

n>0 

6.2 Reception Areas 

Let p(A, r, be the probability to receive a signal at distance r with SIR at least equal to 8. 
Therefore, we have 

p(\,rj,a) = p(w(\)< r jpj 

and the average size of the reception area around an arbitrary transmitter with SIR at least equal 
to 8 is 



a(X,8, a) = 2n J p(X,r, 8,a)rdr . 



Because of the obvious homothetic invariance, we can also write p(A, r, 8, a) = ry/\, 8, a). 
Therefore, a(\,8,a) = ja(l, 8,a). Let a{l,B,a) = 2ir JP(W(1) < L ^)rdr. Assuming that 
w(9, 1) = exp(— CO 1 ) and using integration by parts, we get 



a(l,8,a) = ,aj^p(w(l) = ^) j^dr . 



We use the fact that P(W(1) — x) — ^- J c w(6,l)e ex d6, where C is an integration path in the 
definition domain of w(9, 1), i.e., parallel to the imaginary axis with dt(9) > 0. By changing the 
variable x = (Br")' 1 and inverting integrations, we get 



<7(l,£,a) = ^ fw(9,l) [ e 9x (8x)-^dx 
* l Jc Jo 

= ±p-ir(i- 1 )jw(e,i)(-er- 1 d0. 
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Using w(6, 1) = exp(— CO 1 ), and deforming the integration path to stick to the negative axis, we 
obtain 

0, a) = /T 7 r(l - 7 ) / exp(-C^)^- 1 ^ 

2* Jo 

-, r(i- 7 ) 



= sin(7T7)/3 , 
67 

Using the expression for C, we have 

a(l,/3,a) = /3 7 . 

7T7 

Therefore, the average size of the reception area around an arbitrary transmitter i with SIR at least 
equal to satisfies the identity 

1 sin(-7r) 2 , 
<Ti (\,P, a ) = -—^-J-l3— . (8) 

A —71 

a 



The reception area ai(X,0,a) is inversely proportional to the density of transmitters A and the 
product A<Tj(A, 0, a) is a function of and a. We notice that when a approaches infinity, <7j(A, 0, 00) 
approaches 1/A. This is due to the fact that when a is very large, all nodes other than the closest 
transmitter tend to contribute as a negligible source of interference and consequently the reception 
areas turn to be the Voronoi cells around every transmitter. This holds for all values of 0. The 
average size of Voronoi cell being equal to the inverse density of the transmitters, \, we get the 
asymptotic result. Note that when grows as exp(0(a)), we have <7j(A,/3,a) ~ \ exp(— - log(/3)), 
which suggests that the typical SIR as a — >■ 00 is of the order of exp(0(a)). 

Secondly, as a approaches 2, o"j(A,/3,2) approaches zero because sin(^7r) approaches zero. Indeed, 
the contribution of remote nodes tends to diverge and makes the SIR approach to zero. This 
explains why <7j(A, 0, 2) approaches to zero for any fixed value of 0. 

6.3 Capacity 

In this case, the analytical expressions (jl]) and OH]) lead to 

c(z, 0, a) = E(N(z, 0, A)) = a(l, 0, a) . (9) 



7 Evaluation and Results 



In order to approach an infinite map, we perform numerical simulations in a very large network 
spread over 2D square map with length of each side equal to 10000 meters. 
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7.1 Grid Pattern Based Schemes 



In case of grid pattern schemes, transmitters are spread over this network area in square, hexagonal 
or triangular pattern. For all grid patterns, we set d equal to 25 meters although it will have no 
effect on the validity of our conclusions as capacity, c(z,fl,a), is independent of A. To keep away 
edge effects, we compute the size of the reception area of transmitter i, located in the center of the 
network ) = (0, 0). The network area is large enough so that the reception area of 

transmitter i is close to its reception area in an infinite map. A depends on the type of grid and it is 
computed from the total number of transmitters spreading over the network area of 10000 x 10000 
square meters. In our numerical simulations, we set St = 0.01. 

7.2 Medium Access Schemes Based on Exclusion Rules 

7.2.1 Node Coloring Based Schemes 

Performance of node coloring based schemes is analyzed, via simulations, using the model specified 
in §5.11 We set d equal to 25 meters. 

7.2.2 CSMA Based Schemes 

In order to evaluate the capacity of CSMA based schemes, we perform simulations using the model 
specified in §5.21 The value of carrier sense threshold, 9, is set equal to 1 x 10~ 5 . 

7.2.3 Simulations 

We consider that nodes are uniformly distributed over the network area of 10000 x 10000 square 
meters. Simultaneous transmitters, in each slot, are selected according to the model of each medium 
access scheme. Considering the practical limitations introduced by the bounded network area, we 
use the following Monte Carlo method to evaluate a(\,(3,a). We only compute the size of the 
reception area of a transmitter located nearest to the center of the network area and cr(A, 0, a) 
is the average of results obtained with 10000 samples of node distributions. Similarly, A is also 
the average of the density of simultaneous transmitters obtained with these 10000 samples of node 
distributions. Note that the models of medium access schemes select simultaneous transmitters 
randomly and transmitters are uniformly distributed over the network area. Therefore, using Monte 
Carlo method, i.e., a large number of samples of node distributions and, with each sample, only 
measuring the reception area of a transmitter located nearest to the center of the network area gives 
an accurate approximation of er(A, 0, a) in an infinite map with given values of d or 9. 

It can be argued that, in case of CSMA based scheme, density of simultaneous transmitters is higher 
on the boundaries of the network area, because of lower signal level of interference, as compared to 
the central region. The network area is very large and we observed that the difference, in spatial 
density of simultaneous transmitters, on the boundaries and central region is negligible. We also 
know that the capacity, c(z, (3, a), is independent of A which depends on d or 9. However, if the node 
density is very low, it will also have an impact on the packing (density) of simultaneous transmitters 
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(a) /3 is varying and a is fixed at 4.0. 
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(b) f3 is fixed at 10.0 and a is varying. 



Figure 5: Capacity, c(z,(3,a), of grid pattern (triangular, square and hexagonal), node coloring, 
CSMA and slotted ALOHA based medium access schemes. 

in the network. Therefore, A should be maximized to the point where no additional transmitter can 
be added to the network under given values of d or 9. This can be achieved by keeping the node 
density very high, e.g., we observed that the node density of 1 node per square meter is sufficient 
and further increasing the node density does not increase A. In order to keep away the edge effects, 
values of d or 9 are chosen such that A is sufficiently high and edge effects have minimal effect on 
the central region of the network. 

7.3 Slotted ALOHA Scheme 

In case of slotted ALOHA scheme, capacity, c(z,(3,a), is computed from analytic expressions (|HJ) 
and Q. 
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7.4 Observations 



The values of SIR threshold, 0, and attenuation coefficient, a, depend on the underlying physi- 
cal layer or system parameters and are usually fixed and beyond the control of network/protocol 
designers. However, to give the reader an understanding of the influence of these parameters on 
the capacity, c(z,/3,a), of different medium access schemes, we assume that these parameters are 
variable. 

Figure G^a) shows the comparison of capacity, c(z,(3,a), with grid patterns, node coloring, CSMA 
and slotted ALOHA schemes with j3 varying and a = 4.0. Similarly, Fig. [5(b) shows the comparison 
of these medium access schemes with f3 = 10.0 and a varying. We know that as a approaches infinity, 
reception area around each transmitter turns to be a Voronoi cell with an average size equal to 1/A. 
Therefore, as a approaches infinity, c(z,/3,at) approaches 1. For slotted ALOHA scheme, ([8]) and 
fl9]) also arrive at the same result. For other medium access schemes, we computed c(z,(3,a) with 
a increasing up to 100 and from the results, we can observe that asymptotically, as a approaches 
infinity, c(z, (3, a) approaching 1 is true for all schemes. 

7.4.1 Optimal Local Capacity in Wireless Ad Hoc Networks 

From the results, we can see that the maximum capacity in wireless ad hoc networks can be obtained 
with triangular grid pattern based medium access scheme. In order to quantify the improvement in 
capacity by triangular grid pattern scheme over other schemes, we perform a scaled comparison of 
triangular grid pattern, slotted ALOHA, node coloring and CSMA based schemes which is obtained 
by dividing the capacity, c(z, /3, a), of all these schemes with the capacity, c(z,(3,a), of triangular 
grid pattern scheme. Figure |6] shows the scaled comparison with (3 and a varying. It can be 
observed that triangular grid pattern scheme can achieve, at most, double the capacity of slotted 
ALOHA scheme. However, node coloring and CSMA based medium access schemes can achieve 
almost 85 ~ 90% of the optimal capacity obtained with triangular grid pattern scheme. 

7.4.2 Observations on Node Coloring Scheme 

Triangular grid pattern based medium access scheme can be visualized as an optimal node color- 
ing which ensures that transmitters are exactly at distance d from each other whereas, in case of 
random node coloring, transmitters are selected randomly and only condition is that they must be 
at a distance greater or equal to d from each other. The exclusion region around each transmitter 
is a circular disk of radius d/2 with transmitter at the center. Note that the disks of simultane- 
ous transmitters shall not overlap. The triangular grid pattern can achieve a packing density of 
n/^/l2 « 0.9069. The packing density is defined as the proportion of network area covered by the 
disks of simultaneous transmitters. On the other hand, random packing of disks, which is the case 
in random node coloring, can achieve a packing density in the range of 0.54 ~ 0.56 only [TTjl30"]. 
We have seen in the results that even this sub-optimal packing of simultaneous transmitters by 
random node coloring scheme can achieve almost similar capacity as obtained with optimal packing 
by triangular grid pattern based medium access scheme. 



22 



1 1 


1 1 1 1 1 1 1 

Triangular grid - 




Node coloring " 




CSMA - 




Slotted ALOHA - 










1 1 1 1 1 1 1 1 1 



10 20 30 40 50 60 70 80 90 100 
SIR threshold (P) 

(a) K is varying and a is fixed at 4.0. 



2 

1.75 
1.5 
1.25 

1 

0.75 
0.5 
0.25 




Triangular grid 
Node coloring 
CSMA 
Slotted ALOHA 




10 20 30 40 50 60 70 80 90 100 
Attenuation coefficient (a) 

(b) K is fixed at 10.0 and a is varying. 



Figure 6: Scaled comparison of triangular grid pattern, node coloring, CSMA and slotted ALOHA 
based schemes. 



7.4.3 Observations on CSMA Based Scheme 

We observe that the capacity with CSMA based scheme is slightly lower (by approximately 3%) 
as compared to node coloring based medium access scheme and this is irrespective of the value of 
carrier sense threshold. The reason of slightly lower capacity with CSMA based scheme is that the 
exclusion rule is based on carrier sense threshold, rather than the distance in-between simultaneous 
transmitters. With carrier sense baed exclusion rule, CSMA based scheme may not allow to pack 
more transmitters, in each slot, that would have been possible with node coloring schemes. In other 
words, CSMA may result in a lower packing density of simultaneous transmitters as compared to 
node coloring scheme. This can also be observed by comparing the densities of SSI and SSIk point 
processes in [11] and also explains the slightly lower capacity of CSMA based scheme as compared 
to node coloring scheme. However, as a approaches infinity, A with CSMA based scheme approaches 
the node density and reception area around each transmitter also becomes a Voronoi cell with an 
average size equal to the inverse of node density. In fact, asymptotically, as a approaches infinity, 
capacity, c(z,K,a), approaches 1. 



8 Current Limitations and Future Extensions 

In future, we will extend this work to multi-hop networks. A medium access scheme which achieves 
higher local capacity should also be able to achieve higher end-to-end capacity in multi-hop networks. 
For example, consider that A is normalized across all medium access schemes to 1. Therefore, higher 
local capacity means higher a(l,(3,a) which has an impact on the range of transmission and the 
number of hops required to reach the destination. The analysis to establish exact bounds on end- 
to-end capacity with different medium access schemes in multi-hop networks will be challenging 
as we will have to take into account the impact of routing schemes on capacity as well as various 
parameters like hop length, number of hops and density of simultaneous transmitters which are 
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interrelated. 



The analysis presented here do not take into account fading and shadowing effects. Some results 
with fading are available, e.g., for Poisson distribution of transmitters [7J[9l[T5]. Our analysis, in 
case of slotted ALOHA, can take into account fading by using the results of [20J. Nevertheless, 
analysis of all medium access schemes, discussed here, under the common framework, such as local 
capacity, is lacking. 



9 Conclusions 



We evaluated the performance of wireless ad hoc networks under the framework of local capacity. We 
used analytical tools, based on realistic interference model, to evaluate the performance of slotted 
ALOHA and grid pattern based medium access schemes and we used Monte Carlo simulation to 
evaluate node coloring and CSMA based schemes. Our analysis implies that maximum local capacity 
in wireless ad hoc networks can be achieved with grid pattern based schemes and our results show 
that triangular grid pattern outperforms square and hexagonal grids. Moreover, compared to slotted 
ALOHA, which does not use any significant protocol overhead, triangular grid pattern can only 
increase the capacity by a factor of 2 or less whereas CSMA and node coloring can achieve almost 
similar capacity as the triangular grid pattern based medium access scheme. 

The conclusion of this work is that improvements above ALOHA are limited in performance and 
may be costly in terms of protocol overheads and that CSMA or node coloring can be very good 
candidates. Therefore, attention should be focused on optimizing existing medium access schemes 
and designing efficient routing strategies in case of multi-hop networks. Note that, our results 
are also relevant when nodes move according to an i.i.d. mobility process such that, at any time, 
distribution of nodes in the network is homogeneous. 



.1 Locating the Starting Point z on the Closed Curve Bounding the 
Reception Area 

The SIR Si(z) at point z should be greater or at least equal to 0. We assume that it is equal to 0. 
Therefore, at point z, we have 





Z— Zi\ 




a 


El 


Z-Zj\ 


—a 



Our aim is to find the coordinates of point z which satisfy the above relation. To simplify computa- 
tion of point z on the closed curve, bounding reception area of transmitter i located at Zi = (xi, yi), 
its y coordinate can be fixed such that z = (x,y) = (x,yi). This reduces the above equation to 





X 






El 


z — Z 





Equation ([TO]) is a function of variable x and can be solved using Newton's Method. 



24 



d r» 

(a) Square and Triangular grid based(b) Hexagonal grid pattern based 
schemes. scheme. 

Figure 7: Geometric representation of finding the first approximation of point z = (xo,yi)- 

Remark: Newton's Method: Given a function f(x) and its derivative f'{x), begin with a 
first guess Xq. Provided the function is reasonably well-behaved, a better approximation 
X\ is := x — 4rp-T- The process is repeated until a sufficiently accurate value is reached: 

J \ x ) 

__ fM 



f'(Xn) ' 

From (HUJ) 

We assume that g = r~ a and h = E r 7" where = \x — Xi\ and rj 
Newton's Method requires the derivative of f[x) which is computed as below. 



\Z-Zj\. 



dx 9 " a r? +2 



d \ -\ z — Zj 
—n = > — a ttt 

dx rf +2 



The first derivative of the function f(x) is, 



1 



3^1 
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Newton's Method also requires first approximation of the root, Xq. An approximate value, closer to 
the actual root, can significantly reduce the number of iterations in Newton's Method. 

In all three types of grid networks, the transmitter closest to i, hereafter referred to as j, lies at 
distance d and hence can give the best estimate Xq. For first approximation xq, we can ignore all 
other transmitters in the network. In this case, 



z - Zi 



\z- Zj \ 



where = \z — Zi\ and r\- — \z — zA and we get 



n 

The location of transmitters i and j and point z in the plane form three corners of a triangle with 
angle 9 equal to in case of square and triangular grid and tt/Q radians in case of hexagonal grid 
layout. Figure [7] shows the location of transmitters and Zj, point z and distances r i: rj and d. 
Using above relation between and r, and the Law of Cosines, we get the solution of as 



-B ± VB 2 - AAC 
Ti ~ 2A ' 

where A = 1 — K", B = —2.d. cos(8) and C = d 2 where d is the distance between transmitters i 
and j and is a known parameter of the grid layout. 



Remark: Select positive value ofti as the solution of the above quadratic equation. Using 
Xq = Xi + r i as the first approximate solution in Newton's Method (fTTj) , and after a few 
iterations, we can get a sufficiently accurate value x n+ i which will be the x coordinate 
of the point z. 

The coordinates of point z will be: (x n+ i,yi). 
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